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基于 改进 弦 截 法 的 FastICA 算法 研究 
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摘 要 : 讶 源 分 离 问 题 是 信号 处 理 领域 中 的 研究 热点 之 一 。 众 多 育 源 分 离 算 法 中 国定 点 算法 (FastICA) 因 其 收敛 速度 
快 而 备 受 关注 , 但 是 FastICA 算法 的 收 伊 性 易 受 初始 解 混 矩阵 的 初 值 选择 影响 。 针 对 FastICA 算法 的 不 足 进行 了 改进 ， 
引入 梯度 下 降 法 降低 初 值 选择 敏感 性 ,并且 提出 改进 弦 截 法 , 加快 收 和 你 速度 。 实 验 结 果 显 示 , 基于 改进 弦 截 法 的 FastICA 
算法 与 其 他 FastICA 算法 相 比 ， 不 但 提高 了 算法 的 分 离 性 能 ， 而 且 减 少 了 和 迭代 次 数 ， 增 强 了 收 化 稳 定性 。 所 以 ， 改 进 
的 FastICA 算法 克服 了 初 值 选择 敏感 的 影响 ， 获 得 更 快速 、 更 鲁 棒 的 语音 分 离 性 能 。 
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Research on FastICA algorithm based on improved secant method 
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Abstract: Blind source separation problem(BSS) is one of the hot research in the field of signal processing. In many different 
algorithms for sloving BSS, fixed point algorithm (FastICA) is famous for fast convergence rate. However, its convergence 
performance is easy to affected by initial value selection of the initial demixing matrix. Aiming at the shortcomings of it, this 
paper introduced the gradient descent method to reduce initial value sensitivity, and put forward the improved secant method to 
accelerate the convergence speed. The experimental results show that the improved algorithm compared with the other FastICA 
algorithm, not only improves the separation performance, but also reduces the number of iterations and enhances the convergence 


stability. Therefore, the improved algorithm overcomes the sensitive influence of initial value selection, and achieves faster and 


more robust speech separation performance. 
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盲 源 分 离 (BBS ) 一 直 是 信号 处 理 领 域 中 研究 的 关键 问题 ， 相 比 其 他 ICA 改进 算法 来 说 ， 收 敛 速度 更 快 ， 鲁 棒 性 更 好 四。 
BBS 主要 任务 是 从 多 个 信号 混合 形成 的 观测 信号 中 分 离 出 各 FastICA habc Mp A 中 隐藏 的 独立 
个 不 同 的 源 信号 ,独立 成 分 分 析 ACA) 可 以 有 效 的 解决 这 一 问 ”成 分 非 高 斯 性 最 大 化 的 方法 。 
S, ICA 是 基于 数据 统计 的 信号 处 理 方法 ， 它 能 够 提取 出 隐 FastICA 算法 被 提出 以 后 有 很 多 学 者 对 其 进行 研究 与 改 ; 
藏 在 随机 信号 中 具有 非 高 斯 分 布 和 独立 特性 的 信号 成 分 。ICA ”有 部 分 学 者 针对 FastICA 算法 的 收敛 速度 和 计算 量 进行 优化 。 
在 众多 领域 中 都 取得 了 良好 的 应 用 效果 ， 如 信号 分 离 、 数 字 图 2013 年 Tahir 等 人 0 提出 采用 八 阶 收敛 的 牛顿 法 来 对 FastICA 
像 处 理 、 数 据 分 析 等 方面 中。 算法 进行 优化 ， 很 大 程度 上 提高 算法 的 收敛 速度 ，2017 年 He 
在 20 世纪 80 年 代 早 期 ，Jutten Christian 等 人 在 一 个 神经 ”等 人 0 提出 EFICA 算法 ， 用 两 种 新 的 有 理 数 多 项 式 取代 原来 
生理 学 交流 会 议 上 首次 提出 了 ICA 技术 原型 , 成 功 分离 出 了 两 ”的 非 线性 的 逼近 函数 ， 提 高 了 分 离 性 能 ， 加 快 了 收敛 速度 ， 但 
个 混合 信号 中 。 后 来 许多 改进 的 ICA. 算法 相继 被 提出 ， 如 是 由 于 核心 迭代 公式 仍 是 牛顿 迭代 法 ， 所 以 算法 的 初 值 敏感 性 
Infomax ICA, JADEP!, SOBII9, FastICAU-9I4, 其 中 1997 年 。 并 没有 改善 。 另 外 ， 还 有 部 分 学 者 针对 FastICA 算法 的 收敛 性 
Hyvarinen 等 人 提出 的 固定 点 算法 (Fastpointed ICA, FastICA)， 对 初始 解 混 和 矩阵 的 初 值 选择 比较 敏感 这 一 问题 进行 改进 。2013 
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年 James 等 人 [1 提出 一 种 梯度 迭代 的 ICA 算法 , 使 初 值 矩阵 更 


向 量 w ， 使 得 隐藏 在 混合 信号 中 的 独立 成 分 的 投影 wzz 的 非 高 


快 的 收 化 在 高 斯 噪声 下 取得 更 好 的 稳定 性 ，2014 年 季 策 等 
[31 引 入 超 松 弛 因子 对 初 值 矩 阵 进 行 处 理 , 保 证 收敛 速度 的 前 提 
下 ， 避 免 了 因 初 值 选择 不 当 导 致 的 收敛 不 均衡 现象 ， 2015 年 郭 
伟 等 人 (41 引 入 牛顿 下 山 法 ， 改 善 了 初 值 敏感 性 问题 ， 取 得 了 更 
稳定 的 语音 分 离 效 果 ;2016 年 Ge 等 人 05 提 出 SparseFastICA 算 
法 ， 将 源 信和 号 的 稀 玻 性 作为 FastICA 算法 的 约束 ， 在 磁 共 振 成 
像 (GMRD 中 取得 了 更 好 的 鲁 棒 性 , 但 是 由 于 SparseFastICA 算法 
先 使 用 FastICA 计算 分 离 窍 阵 ， 再 通过 计算 的 源 信号 的 稀 玻 性 
对 分 解 矩 阵 进行 约束 ， 导 致 分 离 效率 不 足 。 经 过 20 多 年 的 发 
展 ， 许 多 改进 的 FastICA 算法 被 提出 ， 但 是 并 没有 哪 一 个 算法 
既 能 提高 收敛 速度 又 能 够 克服 初 值 敏感 性 问题 。 所 以 本 文 尝试 
找到 一 种 既 能 够 更 好 的 降低 初 值 敏感 性 的 影响 又 能 够 加 快 收 敛 
速度 的 方法 。 


1 FastlCA 算法 概述 


FastICA 通过 迭代 不 断 对 解 混 抢 阵 进行 调整 直至 达到 局 部 
最 优 ， 它 是 一 种 快速 寻找 最 优 解 混和 矩阵 的 迭代 算法 中 。FastICA 
算法 主要 包括 两 部 分 ， 数 据 预 处 理 和 独立 成 分 提取 。 

1.1 数据 预 处 理 

数据 预 处 理 通 常 包括 中 心 化 处 理 和 白化 处 理 。 中 心 化 处 理 
是 从 观测 向 量 x 中 减 去 信号 的 平均 值 ， 使 得 观测 向 量 成 为 零 均 
值 向 量 ， 也 称 为 中 心 化 数据 ， 中 心 化 主要 是 为 了 对 混合 信号 进 
行 简化 ， 同 时 降低 噪声 的 影响 。 白 化 处 理 主要 是 对 观测 信和 号 的 
零 均 值 向 量 的 协 方差 矩阵 进行 特征 值 分 解 ， 得 到 中 心 化 数据 的 
化 向 量 z 。 白 化 是 为 了 去 除 中 心 化 数据 的 相关 性 ， 同 时 进行 
降 维 处 理 ， 减 少 计 算 量 ， 一 般 用 主 成 分 分 析 PCA) 进行 白化 
处 理 。 
1.2 独立 成 分 提取 
KURA Xe dE T- ARKE FastICA 算法 ， 算 法 是 以 
PRE Z2 RA BR AE e PE DRE EEUU St, re XOU i s 0] 
Aitte HARAN FATAR: 

JO) ~ pl E(GO) - E(GG)H (1) 
其 中 : JO RREZE, O ERA EEM KN 2808) 
随机 变量 ， 2 表示 一 个 常 系数 ，F{.} 表示 对 数据 求 均值 ，v 是 
具有 零 均 值 和 单位 方差 的 标准 化 的 高 斯 随机 变量 。 为 了 得 到 更 
为 鲁 棒 的 估计 ， 通 常 函数 GC) 选 取 随 自 变量 增长 函数 值 增长 不 
是 太 快 的 函数 ， 被 证 明 非 常 有 用 的 函数 有 00; 


1 
G(y)= zB cosh(a, y) (2) 
Los op ey 
G,(y)= QU 37 (3) 
Go) = 了 (4) 
3C 73? 


通常 a e[L2], a, «1, 
FastICA 算法 可 以 简单 理解 为 找到 一 个 方向 ， 即 一 个 单位 


斯 性 达到 最 大 值 的 迭代 算法 。w'z 的 近似 负 焙 的 极 大 值 通常 在 
E{G(w"z)} 的 极 值 点 处 取得 。 根 据 拉 格 朗 日 条 件 ， EGO 2) 在 
Etw'zy- wl =1 条 件 下 的 极 值 ,是 在 拉 格 朗 日 函数 的 梯度 为 零 
的 点 处 取得 。 所 以 求 负 Jaw) 最 大 值 问 题 可 以 转换 为 求 


E{G(w'z)} 的 极 大 值 问题 ， Ja 00 近似 表示 为 

VeOo)= ElGo2)) « p]w[| -D (5) 
其 中 :B 是 拉 格 朗 日 乘 子 。 对 上 式 两 边 求 偏 导 ， 令 但 导 等 于 零 求 
极 值 点 ， 可 得 目标 函数 Fow) : 

F(w) = E(z g(w"z)} « Bw (6) 


其 中 : g()=G()。 用 牛顿 法 求解 方程 式 (6) ， 寻 找 目标 函数 的 极 
值 点， 牛顿 迭代 法 如 下 所 示 : 


fæ) 
i (7) 


所 以 就 得 到 了 核心 的 迭代 算法 : 


_ , _ Efzg(wz)}+Bw, 
ED (8) 


其 中 : eO 是 8() 的 导数 ，B=Efwizg0wiz)} ， 对 wn 进行 归 一 化 处 
里 ， 提 高 迭代 值 的 稳定 性 。 简 化 后 就 得 到 基于 负 精 的 FastICA 
算法 的 基本 迭代 公式 ; 


(9) 


- E(zg(w1z)) -Elg (wz)}w, 


Wn 
Wan = wanesl 


得 到 最 佳 解 混和 矩阵 w ， 经 过 Y=wx ， 使 得 分 离 矩 阵 
Y «Dy v, RRE i87, 2) isi m y I; BASS T ELA) 
离 信号 ”能 够 尽 可 能 地 逼近 源 信号 。 

由 FastICA 原理 可 知 ， 算 法 依赖 于 初始 解 混 矩 阵 的 选择 和 
迁 代 公式 的 计算 ， 所 以 存在 以 下 两 个 问题 [29: 

a) 初 值 敏感 性 问题 ，FastICA 算法 随机 选取 初始 解 混 和 矩阵 
wo 时， 只 有 we 靠近 收敛 极 值 Y 时 才能 保证 更 好 的 收敛 ， 所 以 
初 值 选 取 会 影响 FastICA 的 收敛 稳定 性 ; 

b) 迭代 过 程 计 算 量 大 的 问题 ， 每 步 欠 代 都 要 计算 gn; 和 
8 (wz) ， 不 论 是 基于 目标 函数 的 求 导 运 算 ， 还 是 雅 可 比 矩 阵 
E(g (w 2)+B 运 算 ， 都 使 得 整个 迭代 过 程 的 计算 量 非 常 大 。 


2 ”改进 弦 截 法 原理 


为 了 减少 迭代 过 程 中 的 计算 量 ， 加 快 欠 代 速度 ， 提 出 一 种 
基于 改进 弦 截 法 的 迭代 方法 ， 避 免 牛顿 迭代 法 的 求 导 运 算 。 
2.4 弦 截 法 

弦 截 法 又 叫 差 商 法 ， 它 是 在 牛顿 迭 代 法 的 基础 上 进行 的 改 
进 ， 用 函数 f(x) 曲线 上 两 点 GfG) 和 Qs fos 构成 的 弦 的 
斜率 近似 的 代替 fo 在 点 %,fG%,) 处 的 切线 斜率 099， 即 : 
Food 


Xn Xn 


fo) (10) 


牛顿 迭代 法 式 (1-7) 和 式 (2-1) 可 推导 出 迭代 公式 为 : 


_ Le x 
mU Fa fL) e 


201 eV 
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这 就 是 弦 截 法 的 迭代 公式 。 
2.2 ”改进 弦 截 法 
在 文献 [17] 中 对 比 了 几 种 不 同 迭 代 方 法 ， 其 中 弦 截 法 的 收 
敛 阶 数 为 1.618, 虽然 略 低 于 牛顿 迭代 法 的 三 阶 收 敛 阶 数 , XE 
速度 稍 显 不 足 ， 但 是 避免 了 牛顿 迭代 法 每 次 兴 代 都 要 求 导 的 耗 
时 计算 .为 了 克服 弦 截 法 收敛 阶 数 不 足 ,迭代 速度 略 慢 的 缺点 ， 
本 文 提出 改进 弦 截 法 。 这 样 既 能 加 快 迭 代 速 度 ， 又 能 减少 每 次 
迭代 过 程 中 的 计算 量 。 改 进 艾 截 法 如 下 : 


Z 2 aaa) 
WT F- Fix 
Cf) (12) 


n4 


2.9 ”改进 弦 截 法 的 收敛 阶 数 

参考 一 阶 弦 截 法 的 收敛 阶 数 的 证 明 0™! 和 ,证 明 改进 弦 截 法 
的 收敛 阶 数 为 2.414， 证 明 过 程 如 下 : 

弦 截 法 是 牛顿 迭代 法 的 变形 ， 根 据 牛 顿 法 的 定义 ， 已 知 方 
程 f(x)=0 有 实 根 ， 设 迭代 序列 {%} 收敛 于 方程 的 根 * ，% 为 每 
次 迭代 的 值 ， 令 e mx -x 表示 第 "次 迭代 的 误差 ， 存 在 实 常数 


p>1 及 kzx0 flfilim ts -k, MERR 


n>% e 


Jj ix) 3& P rli Sk. R 


积 误差 为 es =ke,” , 其 中 ?为 收敛 阶 数 。1 
式 得 


式 (12) 的 第 一 个 分 


= f, )2 1 =f (ka )x, 
mT PO) f Gs) (13) 


结合 迭代 误差 6,= 鸭 -x Hex cx 可 得 


f Ge, — f G, e, 
mT POS) fs) (14) 


根据 拉 格 朗 日 中 值 定理 可 知 ， 在 % 与 x 之 间 存 在 x=5 使 


得 
re- LE 10) (15) 
所 以 有 如 下 变形 公式 : 
f(%)=e,f (&) (16) 
近似 地 有 
Taea, C (17) 
把 式 (16) (17) 代入 式 (15) 可 得 
8, E ef ($)6,, -euf (S,, Je, -e664 f (&) -f (8,4) (18) 
fO) - fo) fO) - fou) 
简化 公式 可 得 
Gn = ke,e, (19) 
其 中 : 大 为 不 等 于 零 的 实数 , 由 定义 知 es =ke, 根据 这 两 个 公 
式 : 
Z =ke? =ke e, = ke, Ê )? kk "e, ” (20) 
所 以 有 如 下 等 式 : 
p=1+ l (21 
» ) 


HX p21, HELA p-(1e5)/2 1618 ， 这 就 是 弦 截 法 的 收敛 
阶 数 。 由 式 〈12) 的 第 二 个 分 式 可 得 ; 
LEa fC g e LESO 
EF) A Faa) FG) 9o) 
所 以 可 以 近似 地 表示 为 
Cnn = KE se, Q3) 
其 中 : 为 非 零 实 数 ， 根 据 式 (20) 可 得 : 
e =k" (Ke,e, e, - k ke,'e, , xu y Ce? (24) 
Hp: C 为 非 零 实 数 ， 由 定义 知 es = 如 ”， 根 据 这 两 个 公式 得 
p-24l (25) 
P 


所 以 p=1+V2=s2414， 改 进 弦 截 法 至 少 
收敛 阶 数 明显 大 于 牛顿 法 的 三 阶 收敛 阶 数 。 


3 ”基于 改进 弦 截 法 的 FastlCA 


3.1 优化 初始 解 混和 矩阵 

为 了 使 初始 解 混和 矩阵 尽快 收敛 ， 减 小 初 值 敏感 性 的 影响 ， 
参考 梯度 法 ?1 ， 引 入 梯度 下 降 法 来 对 初 值 矩 阵 进行 优化 。 梯 
度 下 降 法 又 称 最 速 下降 法 ， 通 常 求 fo 的 极 小 值 点 x ， 一 般 从 
x, 出 发 , 找到 一 个 方向 PA), A xa = 为 +Mp0) ,其 中 为 迭代 步 
长 ， 也 就 是 修正 参数 ， 得 到 梯度 下 降 法 的 目标 函数 : 


有 2.414 阶 收敛 ， 


F&a) = min f x, +Ap(%,)) (26) 
结合 FastICA 算法 式 (6) 可 知 : 
p(w) - Etz gow'z)) (27) 
每 次 迭代 的 修正 参数 可 以 表示 为 
OE{z g(w'z)) 0 
Ow, 
dec 0 m 0 (28) 
" E(z g(w'z)) 


所 以 可 得 w 的 迭代 公式 : 
Win =W, +Àp(w,) =w, +AE{z g(w, z)} (29) 

Fir OARA FE RE RI E IRA 

a) 随机 选择 初始 解 混和 矩阵 w ; 

b) 计算 Elzgowz) 在 % 出 的 负 梯 度 值 ， 计 算 入 s 

c) 更 新 Ww : 

d) WR w 收敛 ， 即 ws -mw|<s ， 独 立成 分 属于 高 斯 分 布 ， 
根据 正 态 分 布 3-s 原则 取 € — 0.00135 ， 则 w,, 即 为 局 部 最 优 的 初 
W, FEER. A, REPER c) 继续 迭代 。 
通过 梯度 下 降 法 降低 初 值 选 择 的 敏感 性 ， 对 随机 初始 矩阵 
进行 修正 ， 能 很 好 的 改善 FastICA 的 初 值 敏感 性 问题 。 另 外 ， 
梯度 下 降 法 的 初期 迭代 速度 很 快 ， 能 加 快 FastICA 算法 的 收敛 
速度 ， 增 强 了 算法 收敛 的 稳定 性 。 
3.2 ”基于 弦 截 法 的 FastICA 

为 了 减少 FastICA 算法 计算 量 ， 文 献 [22] 提 出 


w, =w, +AE{z g(w,'z)) 


NW 


种 基于 牛 
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顿 修 正法 的 FastICA 算法 ， 很 大 程度 上 减少 了 和 迭代 过 程 中 的 计 
算 量 ， 这 给 笔者 带 来 了 启发 。 结 合 目标 函数 式 (6) 和 弦 截 法 式 
d1) ， 可 以 得 到 的 关于 解 混和 矩阵 Y 的 迭代 形式 : 


_ E(w, W, Wna) 
" Fw,)—F(w,) 


Wnn — VW, 


G0) 


所 以 得 到 基于 弦 截 法 的 FastICA 算法 (SFastICA) 的 迭代 公 
X: 


wn = Elzgw,z))w, , - E(zg ws az))w, (1) 


At TUAAREH] FastICA 算法 有 两 个 初始 迭代 矩阵， 和 标准 


分 离 效 果 ， 分 别 对 FastICA、EFICAUN、SparseFastICAL5] (S- 
FICA) 、SFastICA 和 ISFastICA 五 种 算法 进行 实验 比较 。 本 文 
使 用 Two-talker datap32 数据 集 中 三 个 不 同人 的 语音 
(s2 bbbs5p.wav. s5 bbiy7s.wav. s6 bbwz7n.wav) 进行 实验 。 
首先 读 取 三 个 干净 的 语音 信号 ， 然 后 乘 上 一 个 随机 生成 的 混合 
系数 和 矩阵， 得 到 混合 信号 ， 之 后 分 别 用 五 种 算法 对 混合 信号 进 
行 分 离 ， 分 离 得 到 三 个 估计 信号 ， 最 后 分 别 统计 五 种 算法 的 分 
离 信 号 的 信 噪 比 、 和 迭代 次 数 、 分 离 时 间 ， 并 分 析 算 法 的 收敛 特 
性 和 对 初 值 选择 的 敏感 程度 。 


FastICA 算法 一 样 易 受 初 值 选择 影响 。 但 是 引入 梯度 下 降 法 对 
初 值 矩阵 进行 优化 ， 将 弥补 SFastICA 算法 收敛 速度 不 足 的 缺 
点 。 所 以 ，SFastICA 算法 既 能 减少 初 值 敏感 性 又 能 在 一 定 程度 
上 提高 收敛 速度 和 稳定 性 。 
3.8 ”基于 改进 弦 截 法 的 FastICA 

收敛 阶 数 越 高 从 代 收敛 速度 越 快 ， 为 进一步 提高 弦 截 法 的 
收敛 速度 ， 本 文 提出 基于 改进 弦 截 法 的 FastICA 算法 
(ISFastICA) 。 改 进 弦 截 法 在 2.3 Bou Sx Er C 
2.414， 具 有 比 牛 顿 迭代 法 更 高 的 收敛 阶 数 ， 收 敛 速度 明显 高 于 
牛顿 兴 代 法 。 结 合 目 标 函 数 式 (6) 和 改进 弦 截 法 式 (120 , "T 


得 到 迭代 公式 ， 
_ (w, -w, Fn) 
nd 7 W,— Fw )— Fr, ) 
n n-l 
(Ws - WF) 682) 
Wa 7W,— E 
FO, D- F(v,) 


所 以 可 以 推导 出 基于 改进 弦 截 法 的 FastICA SEXES INA: 


WwW, = BlizgOw; z))w, , —E{zg wz) }w, 
wn 7 E{zg (Wa)  z))w, — E{zg Ov; z) Ww, (33) 


Wia” wall". | 


所 以 基于 改进 弦 截 法 的 FastICA 算法 步骤 如 下 : 


`; mosh D INS 
a) 对 观测 信号 进行 中 心 化 处 理 : xx 251 


b) 对 中 心 化 数据 进行 白化 ， 
z =Vx = ED"E"x ; 

c) 随机 选择 两 个 具有 单位 范 数 的 初始 化 向 量 w。 和 ww : 
Inl - ef =1， 优 化 初 值 矩阵 ; 

d) 更 新 w : 将 优化 的 解 混和 矩阵 代入 式 (30 中 得 到 w ; 

e) WR w dk, 即 wa-w,|<8 ,根据 正 态 分 布 4-o 原 则 取 
z=0.0000315 ， 则 算法 收 化 ， 停 止 迭 代 。 否 则 返回 步 又 qd); 

f) 将 收敛 后 的 最 优 解 混 矩阵 ” 与 观测 信号 x 相 乘 ， 得 到 源 
信号 的 估计 信号 了 s 

ISFastICA 算法 从 收敛 特性 上 明显 比 SFastICA 算法 和 标准 
FastICA 算法 具有 更 大 的 优势 。 它 不 但 能 降低 初 值 敏感 性 影响 ， 
而 且 收 敛 阶 数 更 高 ， 收 敛 速度 更 快 。 所 以 ISFastICA 算法 的 分 
离 特 性 和 稳定 性 更 好 。 


得 到 白化 向 量 : 


4 ”实验 结果 及 分 析 


为 了 评估 基于 改进 弦 截 法 的 FastICA 算法 (ISFastICA) 的 


4.1 分 离 性 能 对 比 实验 


分 离 实验 用 信 噪 比 作为 衡量 分 离 效果 的 好 坏 标准 ， 信 噪 比 
可 以 简单 描述 为 估计 信号 中 的 有 效 成 分 和 噪声 成 分 的 比值 ， 一 


般 是 信 噪 比 大 的 说 明 分 离 效果 更 好 2。 信 噪 比 定 义 如 下 : 
MU 
SNR -10*1og,, | 一 过 


Ys 
i-0 


(4) 


其 中 :表示 干净 的 源 信号 ，$ 表示 分 离 的 估计 信和 号。 本 节 进 行 
100 次 分 离 实验 , 对 五 种 算法 的 分 离 实验 的 信 噪 比 、 和 欠 代 次 数 、 
分 离 时 间 进 行 统计 ， 求 均值 ， 如 表 d 所 示 。 

del 实验 参数 对 比 


算法 信 噪 比 /dB RAA 。 分 离 时 间 /s 
FastICA 40.5719 18 0.9721 
EFICA 42.4302 16.4 0.8724 
S-FICA 41.8616 14.6 0.9876 
SFastICA 41.6004 16 0.8803 
ISFastICA 42.1221 n 0.7318 


表 1 可 知 , 对 比分 离 信号 的 信 噪 比 , ISFastICA 算法 的 分 
离 效果 优 于 S-FICA 算法 优 于 SFastICA 算法 优 于 FastICA 算法 ， 
略 输 于 EFICA 算法 。ISFastICA 算法 的 信 噪 比比 FastICA 算法 
提高 了 3.82%, SFastICA 算法 的 信 噪 比比 FastICA 算法 提高 了 
2.54%. XE TQ, ISFastICA 算法 少 于 S-FICA 算法 少 于 
SFastICA 算法 少 于 少 于 EFICA 算法 少 于 FastICA 算法 。 
ISFastICA 算法 迭代 次 数 比 FastICA 算法 降低 了 38.8996, 
SFastICA 算法 迭代 次 数 比 FastICA 算法 降低 了 11.1196. 对 比分 
离 时 间 , ISFastICA 算法 快 于 EFICA 算法 快 于 SFastICA 算法 快 
于 FastICA 算法 快 于 S-FICA 算法 。ISFastICA 算法 分 离 时 间 比 
FastICA 算法 提高 了 24.72%,SFastICA 算法 分 离 时 间 比 FastICA 
算法 提高 了 9.44%。 所 以 ISFastICA 算法 不 论 是 分 离 效果 还 是 
分 离 效率 都 有 很 大 的 提升 。 
4.2 算法 收敛 性 对 比 实验 

读 入 三 个 人 的 源 信号 ， 图 1 是 的 
， 横 轴 为 采样 点 ， 纵 轴 为 波形 幅 值 。 

对 三 源 信号 进行 线性 混合 ， 乘 以 一 个 随机 生成 混合 系数 矩 
阵 ， 得 到 三 个 混合 信号 的 波形 图 ， 如 图 2 所 示 。 
使 用 五 种 算法 分 别 对 混合 信号 进行 分 离 估计 ， 得 到 不 同 算 
法 的 分 离 结果 ， 如 图 3 所 示 。 


净 语 音信 号 的 采样 波形 


NS 


加 
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答 入 声音 1 Rari 
05 T T T 05, T T T T T 
Roe Ae cede m m 二 
Be 1 2 : r 7 1 1 1 1 1 
se j 15 2 2 E zd M hi H M 05 1 15 2 25 3 35 4 45 5 
Mraz : 混合 声音 2 xi 
5 05 T T T T 05 T T T 
45 L T : " t 1 1 1 1 L 
d e , 13 E e 3 E * A: : 0 05 1 15 2 25 3 35 4 45 5 
输入 声音 3 hi irf xv 
05 - r T [2 : r - 
a5 L 1 L 1 .05 L n n L 1 
0 05 1 15 2 25 3 35 4 45 5 9 os 1 15 2 25 3 35 4 45 M 
f 10 
采样 点 bd 采样 点 3 
ZA [zi Y VB. vr. [22 EL LN 局 yH AJE 422 ET 
图 1 个 不 同人 的 纯净 语音 信号 图 2 三 个 不 同人 的 混合 语音 信号 
FastcA-1 EFICA-1 
10 T T T T T 10 T T 
- mM 
o Ars lien tmm i omm of- Aide — 
" I : 2 : 3 i i 95 1 15 2 25 3 35 4 45 5 
o 05 t 15 2 25 3 35 4 45 5 S n E ie 
10 x 10 
FastiCA 2 EFICA-2 
10. 10 
E E 
X oF- 7 MAR M mime x 0 " 
a Š T 
E 
"s 1 n 1 n -10 1 
[] 05 1 15 2 25 3 35 4 45 5 05 1 15 2 25 3 35 4 45 5 
10 x 10 
o Fastca-3 EFICA-3 
10 
ET 1 | - ame e 
| j 
"S 05 1 15 2 25 3 35 4 45 5 -10 k 
村 0.5 1 15 2 25 3 35 4 45 5 
x 
采样 点 采样 点 x 10* 
sd du NES 
(a) FastICA 算法 效果 (b) EFICA 算法 效果 
S-FICA-1 
ag : Sferen P — 
j 
Hmm — I — " 
9 Afi rile t o ov 中 ———— 
E 1 i 1 1 1 1 
0 os 1 15 2 25 3 35 4 45 5 ^ P" 
vw $——39——1—— 359 —1- —5-——1 —4—1—$—3 os EE ss 
x h 
" S-FICA-2 PETS id ISFastCA-2 nici 
o 
" : a 
HE Alti ullo te — m mm £o dfe idi Mijas ina HE jain delle «t laedi tmi 
Fio La ah L s M" 
0 05 1 15 2 25 3 35 4 45 5 tu 1 315 — E —18 —3 —35 —4 —45 —3 9 0205 1 1 23 35 3 3 4 4 5$ 
n ; sranna " rancas 
SFICA-3 1 
10 
a " : E n " " 
5— 98 —1 18 —3 38 3 38 ^4 ^48 ^s $25 worth od 5 3 od 8539 358—158 3658384 
un "m 采样 点 采样 点 


(c) S-FICA 算法 效果 (d) SFastICA 算法 效果 (e) ISFastICA 算法 效果 
图 3 五 种 算法 分 离 结果 


统计 五 种 算法 分 离 时 的 解 混 矩阵 迭代 差 值 和 迭代 次 数 之 间 ”的 初 值 矩阵 ， 然 后 统计 分 析 不 同 算法 的 迭代 次 数 变化 ， 如 表 2 
的 关系 ， 评 价 算法 的 收敛 特性 ， 和 迭代 差 值 下 降 的 越 快 ， 则 算法 MR 
收敛 速度 越 快 。 收 敛 特 性 的 三 维 立 体 柱状 统计 图 如 图 4 所 示 ， 


从 外 向 内 分 别 代表 FastICA、EFICA、S-FICA、SFastICA 和 le 


ISFastICA 算法 , 用 由 浅 到 深 色 调 的 方形 柱 体 表 示 , 纵 轴 用 柱 体 053 A 

(Es REESE IARE, Bp II OC. zd I 
由 图 4 可 以 看 出 , FastICA、EFICA、S-FICA、ISFastICA、 - 
SFastICA、FastICA 迭代 次 数 分 别 为 18、13、14、16、11。 明 显 i 
0 


可 以 看 出 ， 和 迭代 收敛 速度 由 快 到 慢 分 别 是 ISFastICA, EFICA, 

S-FICA、SFastICA、FastICA 算法 。 所 以 ，ISFastICA 收敛 性 能 

最 好 。 

4.3” 初 值 敏感 性 对 比 实验 送 代 次 数 S161718 
为 了 对 比 初始 解 混和 矩阵 的 初 值 选择 对 算法 收敛 性 的 影响 ， 

对 三 个 人 的 混合 信号 进行 10 次 分 离 实验 ， 每 次 随机 产生 不 同 


图 4 算法 收敛 性 对 比 
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RARA 
No. 

FastICA EFICA S-FICA SFastICA ISFastICA 

1 18 16 17 16 11 
2 24 14 14 13 7 
3 23 10 15 15 10 
4 21 13 14 13 8 
5 16 12 13 17 11 
6 20 11 11 15 7 
7 18 14 12 11 9 
8 25 16 13 14 7 
9 22 9 15 13 10 
10 16 11 14 17 7 
均值 20.3 12.6 13.8 14.4 8.7 
方差 10.46 5.24 2.56 3.82 2.90 


EFICA 小 于 


迭代 次 数 统 计 表 2 可 知 ，ISFastICA 迭代 次 数 均值 小 于 
F S-FICA 小 于 SFastICA 小 于 FastICA, 
迭代 次 数 方差 小 于 SFastICA 小 了 


F EFICA 小 于 FastICA， 稍 大 


于 S-FICA。 所 以 ，ISFastICA 算法 迭代 次 数 最 少 ， 


值 解 混 


较 小 。 


为 了 更 直观 地 看 出 每 次 实验 时 ， 随 机 生成 的 初 


代 收 敛 性 能 的 影响 ， 把 每 次 实验 的 迭代 次 数 绘 


区 


5 看 出 ， 每 一 组 ISFastICA 


迭代 次 数 


且 各 组 之 间 的 
解 混 矩阵 ISFastICA 算法 日 


ERI 


直方 


高 低 差 别 不 大 ， 


f 法 选 代 次 数 对 


的 迭代 次 数 直 
反映 出 对 于 不 同 的 初 值 
的 收敛 性 能 最 稳定 ， 也 就 是 说 
ISFastICA 算法 相 比较 其 他 四 种 算法 对 初 值 敏感 性 影响 最 小 ， 

而 且 迭 代 收 敛 速度 最 快 。 


i 


ISFastICA 


对 


f 


HRES E 


4.4 实验 结果 分 析 


对 实验 结果 进行 统计 分 析 ， 对 比 五 种 算法 的 实验 结果 ， 可 


以 看 出 可 以 得 到 如 下 结论 : 


a) 


d 1 可 知 ,基于 改进 弦 截 法 的 FastICA 算法 (ISFastICA) 


分 离 性 能 优 于 EFICA 优 于 S-FICA 优 于 SFastICA 优 于 FastICA 
算法 。 


b) 对 比 算法 收敛 特性 ， 如 图 


4 所 示 , 可 以 看 出 : ISFastICA 


算法 收敛 速度 最 快 。 


种 算法 ， 不 仅 收 敛 速 度 更 快 而 且 收敛 稳定 性 更 好 。 


c) 初 值 敏感 性 对 比 实验 可 知 ，ISFastICA 算法 对 比 其 1 


(E=; 


ISFastICA 算法 对 比 标准 FastICA 算法 , 不 但 降低 了 初 值 选 


m 


择 敏感 性 , 而 且 收 敛 速度 更 快 ， 收敛 稳定 性 更 好 。ISFastICA 算 


法 对 比 SFastICA 算法 , 提高 了 算法 的 分 离 性 能 , 加 快 了 收敛 速 


度 。ISFastICA 算法 对 比 EFICA 算法 虽然 


SEHE (SNR) KA 


不 足 , 但 是 收敛 速度 更 快 、 鲁 棒 性 更 好 。ISFastICA 算法 对 比 S- 


FICA 算法 虽然 收敛 稳定 性 略 
有 了 明显 的 提高 。 


9 不足 ， 但 是 收敛 速度 和 SNR. 都 
对 比 FastICA、SFastICA 和 ISFastICA 三 种 


A 


算法 可 以 看 出 ， 梯 度 下 降 法 的 引入 有 效 地 降低 了 初 值 选择 敏感 


性 的 影响 ， 使 算法 的 迭代 次 数 方差 更 小 ， 


收敛 更 稳定 。 所 以 改 


FREY 


JERE RISE E PE LEERI, SEAN EIA TELE USE RE] E 


矩阵 对 和 迭 
, 如 


图 4-5 所 示 ， 每 一 组 从 左 到 右 分 别 表 示 ISFastICA、EFICA、S- 


FICA, SFastlCA, FastlCA 算法 ， 颜 色 从 浅 到 深 ， 纵 贡 
图 的 高 度 表 示 友 代 次 数 ， 横 轴 表 示 实 验 序号 。 


lI 


直方 


方 


氏 ， 


图 最 


口 FastICA Œ EFICA Œ S-FICA Bi SFastICA Bi ISFastICA 


进 弱 截 法 的 提出 在 提高 分 离 性 能 的 前 提 下 ， 既 能 够 降低 初 值 敏 


感性 ， 又 能 够 加 快 收敛 速度 ， 增 强 收敛 稳定 性 。 
5 HRA 


提出 改进 弦 截 法 ,提高 算法 的 收敛 速度 ,增强 算法 收敛 稳定 性 。 


本 文 引入 梯度 下 降 法 克服 FastICA 对 初 值 选 择 敏感 的 缺陷 ， 


通过 理论 证 明和 实验 验证 ，ISFastICA 算法 收敛 速度 快 于 
FastICA 算法 快 于 SFastICA 算法 。ISFastICA 算法 收敛 稳定 性 
优 于 SFastICA 算法 优 于 FastICA 算法 。 另 外 ,通过 与 当前 最 新 
的 FastICA 的 改 ; 
ISFastICA 算法 的 
提高 。 所 以 ， 基 于 改进 弦 截 法 的 FastICA 算法 不 
度 ， 
算 ,所 以 


is 


| 


算法 EFICA 和 SparseFastICA(S-FICA) 相 比 ， 
性 敛 稳定 性 和 分 离 性 能 分 别 也 有 一 定 程度 的 
晶 提 高 收敛 速 
曾 强 收敛 稳定 性 。 但 是 由 于 改进 的 算法 无 先 验 知识 计 
不 能 够 区 分 和 提取 混合 语音 中 某 一 个 特定 人 的 语音 ， 


而 且 


步 将 进行 深度 聚 类 的 语音 分 离 研究 。 
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